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ON  THE  STABILITY  OF  A  TIME 
DEPENDENT  BOUNDARY  LAYER. 


S.  R.  Otto* 

ICASE,  Mailstop  132c,  NASA  Langley  RC, 
Hampton,  VA  23681-0001,  USA. 


ABSTRACT 


The  aim  of  this  article  is  to  determine  the  stability  characteristics  of  a  Rayleigh 
layer,  which  is  known  to  occur  when  the  fluid  above  a  flat  plate  has  a  velocity  im¬ 
parted  to  it  (parallel  to  the  plate).  This  situation  is  intrinsically  unsteady,  however 
as  a  first  approximation  we  consider  the  instantaneous  stability  of  the  flow.  The 
Orr-Sommerfeld  equation  is  found  to  govern  fixed  downstreeim  wavelength  linear  per¬ 
turbations  to  the  beisic  flow  profile.  By  the  solution  of  this  equation  we  can  determine 
the  Reynolds  numbers  at  which  the  flow  is  neutrally  stable;  this  quasi-steady  approach 
is  only  formally  applicable  for  infinite  Reynolds  numbers.  We  shall  consider  the  large 
Reynolds  number  limit  of  the  original  problem  and  use  a  three  deck  mentality  to  de¬ 
termine  the  form  of  the  modes.  The  results  of  the  two  calculations  are  compared,  and 
the  linear  large  Reynolds  number  analysis  is  extended  to  consider  the  effect  of  weak 
nonlinearity  in  order  to  determine  whether  the  system  is  sub  or  super  critical. 
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No.  NASl-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681 
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§1  Introduction 

j  The  main  thrust  of  boundary  layer  stability  calculations  have  involved  spatiallv 

I 

evolving  problems.  Relatively  little  work  ha.s  gone  into  the  stability  of  temporally 
developing  layers,  one  such  situation  is  the  topic  of  this  paper,  that  of  an  impulsive 
flow  over  a  flat  plate. 

The  calculation  to  determine  the  linear  stability  of  the  system  considered  here  is 
initially  attempted  using  a  quasi-steady  approach  where  the  flow  is  frozen  at  a  certain 
time.  This  yields  an  Orr-Sommerfeld  equation,  the  solutions  of  which  were  discussed 
for  a  variety  of  flows  by  Tollmien  (1929)  and  Schlichting  (1933),  using  homogeneous 
assumptions  in  the  temporal  and  downstream  coordinates.  Unfortunately  the  true 
system  governing  what  we  shall  refer  to  as  Tollmien-Schlichting  waves  is  not  uniform 
for  a  general  boundary  layer.  An  attempt  to  remedy  this  deficiency  was  propo.sed  for 
the  Taylor-Gortler  problem  by  Smith  (1955).  These  ideas  were  extended  by  Bouthier 
(1973)  and  Gaster  (1974),  where  extra  terms  were  added  to  appreciate  for  the  boundary 
layer’s  growth;  although  it  was  noted  by  Gaster  that  this  only  amounted  to  a  successive 
approximation  technique.  It  was  in  the  work  of  Smith  (1979  a)  that  a  more  formal 
asymptotic  treatment  of  the  problem  was  proffered  for  large  Reynolds  numbers.  Smith 
used  a  triple  deck  structure  to  describe  the  linear  development  of  fixed  frequency 
Tollmien-Schlichting  waves  within  a  growing  boundary  layer,  specifically  a  Blasius 
layer.  It  was  found  that  the  boundary  layer  growth  altered  the  higher  order  terms  in 
the  expansions  of  various  wave  characteristics  within  this  structure.  The  inclusion  of 
these  terms  may  be  stabilising  or  destabilising  depending  on  which  quantity  is  taken 
to  be  representative  of  the  perturbations  intensity. 

Some  studies  have  been  made  of  the  stability  characteristics  of  unsteady  bound¬ 
ary  layers,  although  most  of  the  cases  have  involved  temporal  periodicity.  In  these 
problems  either  a  quasi-steady  approach  or  Floquet  theory  can  be  used.  In  Kerczek 
&  Davis  (1974)  the  linear  stability  of  a  Stokes  layer  on  a  flat  plate  is  discussed,  in  this 
problem,  waves  that  develop  quickly  compared  to  the  temporal  scale  of  the  basic  state 
are  considered,  as  is  the  case  in  Cowley  (1986).  In  Hall  (1975)  a  study  was  made  of 

I 

unsteady  flows  about  cylinders,  and  in  Seminara  &  Hall  (1976)  the  Taylor  problem  is 
considered  where  the  inner  cylinder  is  taken  to  oscillate  periodically  with  time.  The 
work  of  Di  Prima  &  Stuart  (1975)  regards  the  Taylor  problem  for  eccentric  cylinders 

I  (which  corresponds  to  a  journal  bearing).  In  Otto  (1993)  ^he  stability  of  a  boundary 
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layer  at  an  infinite  cylinder  is  discussed.  The  cylinder  starts  to  spin  up  and  vortices 
are  generated  by  an  axisymiuetric  element  of  wall  roughness;  these  may  lead  to  some 
kind  of  transition.  In  Duck  (1990)  the  unsteady  triple  deck  equations  were  solved  for 
the  flow  over  a  small  hump  on  a  flat  plate. 

in  the  current  problem  the  basic  state  depentls  on  the  normal  and  temporal  coor¬ 
dinates,  so  it  is  natural  to  solve  an  unsteady  problem.  This  situation  is  fully  parallel 
so  there  is  no  need  to  make  any  assumptions  about  the  spatial  evolution  of  the  basic 
state.  After  a  time  toi  we  may  define  the  boundary  layer  Reynolds  number  as, 

D  Uoail'to)"^  Uootl 
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based  on  the  viscous  layer’s  thickness.  The  velocity  imparted  to  the  fluid  is  f’oo« 
and  its  kinematic  viscosity  is  2v.  In  the  viscous  layer  (which  may  be  shown  to  have 
thickness  (/^t)  ^ ),  the  basic  state  is  found  to  be  similar  in  time  and  is  given  by  an  error 
function.  This  is  accreditable  to  the  large  Reynolds  number  assumption;  at  finite 
Reynolds  numbers  the  unsteady  equations  need  to  be  solved.  We  wish  to  consider  the 
susceptibility  of  this  profile  to  Orr-Sommerfeld  modes  that  evolve  on  short  temporal 
and  spatial  scales.  As  mentioned  previously  the  approximation  is  made  that  we  shall 
‘freeze’  the  basic  profile.  In  an  infinite  Reynolds  number  flow  this  is  formally  correct 
as  there  is  no  way  information  can  propagate;  but  in  the  Orr-Sommerfeld  method 
the  Reynolds  number  is  retained  as  a  finite  parameter.  In  this  paper  we  do  not 
wish  to  discuss  the  legitimacy  of  the  Orr-Sommerfeld  approach,  it  is  just  used  eis  a 
first  approximation  to  the  solution  of  the  linear  stability  problem  at  finite  Reynolds 
number.  For  each  time  to  beyond  a  certain  point  (this  corresponds  to  Reynolds  number 
Reo)^  this  problem  may  be  solved  to  obtain  a  real  wavenumber  and  phase  speed  such 
that  the  wave  is  neutral.  Above  this  critical  Reynolds  number  there  are  two  neutral 
values,  and  as  the  Reynolds  number  increases  still  further  the  behaviour  of  the  wave 
having  the  smaller  wavelength  is  governed  by  a  five  deck  structure,  the  critical  layer 
becomes  distinct  from  the  wall  layer,  due  to  the  size  of  the  phase  speed.  Our  interest 
is  confined  to  the  lower  branch  where  the  structure  may  be  described  using  three 
decks.  An  ample  description  of  this  structure  can  be  found  in  Stewartson  (1974).  The 
bulk  of  the  boundary  layer  is  referred  to  as  the  main  deck,  which  at  leading  order  is 
inviscid  and  linear.  In  the  proximity  of  the  wall  there  is  a  fully  viscous  layer,  this 
region  contains  the  critical  layer  of  classical  Orr-Sommerfeld  theory.  Finally  there 
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is  an  upper  deck  which  is  external  to  the  conventional  boundary  layer  and  enables 
matching  with  the  outer  potential  flow  solution.  A  multiple  scales  technique  is  used 
to  resolve  the  growth  of  the  boundary  layer  with  time. 

The  temporal  growth  of  the  layer  produces  corrections  to  the  eigenproblems  for 
the  wavenumbers  and  phase  speeds  at  a  higher  order.  The  triple  deck  solutions  may 
be  compared  with  the  Orr-Sommerfeld  calculation  in  the  regime  of  common  validity, 
that  is  where  i?eo  >>  1-  In  fact  the  corrections  caused  by  the  growth  of  the  layer  to 
the  ‘steady’  triple  deck  scenario  occur  in  the  third  order  equations,  one  order  earlier 
than  those  in  Smith  (1979  a).  This  is  due  to  the  different  temporal  and  spatial  scales 
over  which  Tollmien-Schlichting  waves  develop. 

The  results  of  this  calculation  may  be  viewed  in  another  sense.  For  a  perturbation 
wavenumber  o,  at  each  time  to  we  have  a  frequency  Q.  When  the  imaginary  part  of  Q 
becomes  zero,  the  wave  is  instantaneously  neutral.  In  the  neighbourhood  of  this  neutral 
time  we  may  develop  a  weakly  nonlinear  theory.  A  Stuart-Watson  amplitude  equation 
can  be  derived,  and  using  this  we  may  determine  whether  the  perturbations  are  sub 
or  super-critical,  that  is  whether  the  nonlinearity  has  a  destabilising  or  stabilising 
effect  on  the  situation.  We  shall  also  consider  the  interplay  between  the  unsteady  and 
nonlinear  effects,  as  was  discussed  for  the  Blasius  problem  by  Smith  (1979  b)  and  Hall 
&  Smith  (1984). 

The  procedure  adopted  in  the  remainder  of  this  paper  is  as  follows,  in  section  2 
the  basic  flow  is  derived  and  the  queisi-steady  problem  is  formulated.  In  section  3  the 
linear  triple  deck  problem  is  formulated  and  solved.  In  section  4  the  weakly  nonlinear 
problem  is  derived  and  the  nonlinear  amplitude  equation  is  given.  Section  5  contains 
details  of  the  numerical  solution  of  the  linear  problem.  Section  6  includes  discussion 
of  the  effect  of  considering  certain  small  but  finite  disturbances,  and  in  section  7  some 
brief  conclusions  are  drawn. 


§2  Derivation  of  the  basic  flow  and  the  formulation  of  the  quasi-steady  problem 


The  problem  is  non-dimensionalised  in  the  usual  way  and  the  unsteady  three  di¬ 
mensional  Navier-Stokes  equations  in  cartesian  coordinates  become  (with  the  temporal 
scaling  t  ~  L/Uoo  for  the  Orr-Sommerfeld  modes,  where  L  is  a  typical  lengthscale,  and 
t  ju  for  the  basic  state). 


du  du  du  du  dp  1 


u. 
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dv  dv  dv  dv 

= 

dw  dw  dw  dw 

■S'+“87  +  ''^  +  "’87  = 

du  dv  dw 
dx  ^  dy  ^  dz 


dp  1  _2 
dy  Re 

dp  1  2 

’ -51  +  57^ 
=  0, 


where  Re  is  the  Reynolds  number  defined  as  UooL/v,  and  is  the  Laplacian  operator 
in  cartesian  coordinates.  The  choice  of  the  lengthscale  L  is  discussed  in  section  5,  but 
it  shall  be  taken  to  be  the  displacement  thickness  of  the  layer,  as  in  Jordinson  (1970). 
In  this  section  disturbances  shall  be  considered  to  be  two  dimensional  to  facilitate 


comparison  between  this  work  and  that  of  Smith  (1379  a),  to  this  end  it  is  assumed 
that  d/dz  =  0  and  u;  =  0. 

We  require  that  the  unsteady  terms  and  viscous  terms  balance  in  the  viscous  layer, 
this  implies  that  this  layer  is  of  thickness  ■  In  the  non-dimensional  coordinates 
this  corresponds  to  Thus  we  introduce  a  boundary  layer  similarity  variable  r/  giveii 

by 


y 

"  =  5T 


(2.1) 


There  is  no  flow  in  the  normal  direction,  as  can  be  seen  from  the  equation  of  continuity 
at  leading  order.  Hence  the  basic  state  is  governed  by 


d'^UB 

5r/2 


+  2t7 


duB 

dt] 


=  0, 


along  with  the  boundary  conditions 


ub  =  0  at  t/  =  0,  and  ub  I  as  7/  — >  00. 
The  solution  of  his  system  is  given  by 


1 


0 


this  is  an  error  function. 

Now  in  the  context  of  an  Orr-Sommerfeld  equation  we  consider  the  instant  sta- 
bility  of  this  profile,  hence  (2.1)  now  becomes  7}  =  y/2tQ  ,  where  to  is  a  constant.  We 
suppose  that  any  perturbation  to  the  flow  is  independent  of  the  order  one  temporal 
scale. 


4 


The  perturbations  are  considered  to  be  travelling  waves  which  evolve  on  the  scale 
L,  thus  the  disturbance  is  proportional  to  exp(z^(x  —  ct)).  The  disturbance  is  taken 
to  be  infinitesimal  so  that  the  resulting  analysis  is  linear.  As  the  perturbation  is  two 
dimensional,  a  streamfunction  formulation  can  be  exploited,  such  that 


dy 


V  =  - 


drp 

dx 


It  should  be  noted  that  terms  proportional  dtl^/dt  =  0  have  been  ignored,  so  that  the 
equation  may  be  solved  as  an  ordinary  differential  equation.  The  governing  equation 
is  the  conventional  Orr-Sommerfeld  equation  namely. 


( 


/  / 


idRe  \dy 


^  -20^  ® 


C)y2 


(2.2) 


The  boundary  conditions  imposed  are  those  of  no-slip  at  the  wall,  and  the  decay  of 
the  disturbance  as  y  ^  oo,  this  implies  the  wave  is  confined  to  the  boundary  layer. 
This  eigenvalue  problem  can  be  solved  numerically  employing  the  method  discussed  in 
Malik,  Chuang,&  Hussaini  (1982).  The  method  uses  a  two-point  fourth  order  compact 
finite  difference  scheme  based  on  a  Euler- Maclaurin  method.  A  stretched  grid  was  used 
so  that  resolution  could  be  retained  without  using  a  prohibitive  number  of  points.  The 
c.-'l  julations  were  checked  by  finding  the  eigenvalues  of  the  adjoint  system  to  (2.2).  The 
critical  values  were  found  to  be 

Rec  =  2968.4,  Uc  =  0.13340,  9c  =  0.4255,  c  =  0.3135, 

where  Q  =  9c.  This  point  represents  the  smallest  value  of  the  boundary  layer  Reynolds 
number  at  which  the  flow  situation  is  linearly  unstable  to  travelling  waves  of  the  type 
considered  herein.  Plots  of  9  and  c  against  Re  are  given  in  figures  1  and  2  respectively. 
These  calculation  were  affected  for  to  =  1/4,  so  that  y  =  rj  at  that  instance. 


Figure  l:Streajnwise  wavenumber  versus  i?e,  Orr-Sommerfeld  calculation 


Figure  2:Phase  speed  versus  Re,  Orr-Sommerfeld  calculation 

§3  Formulation  and  solution  of  the  triple  deck  problem  in  the  linear  regime 

In  this  section  the  structure  of  the  the  linear  disturbances  is  described  using  triple 

_  I 

deck  theory.  As  mentioned  previously  i?e  is  assumed  to  be  large;  and  e  =  i?e  *  ‘C  1  is 
introduced.  In  classical  Orr-Sommerfeld  theory  the  typical  wavelength  of  the  neutrally 
stable  modes  increases  proportionally  to  Re  as  i2e  1,  see  Stuart  (1963).  Stuart 
confirms  this  asymptotic  property  with  a  Reynolds  number  based  on  boundary  layer 
thickness.  Three  distinct  regions  are  considered,  namely  the  upper,  main  and  lower 
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decks,  scaled  normally  on  e^,  e'*,  and  e®  respectively.  Further  temporal  and  spatial 
variables  are  introduced,  that  is  T  and  X,  given  by,  T  =  and  .Y  =  The 

solutions  are  considered  .n  the  three  decks  in  terms  of  expansions  involving  powers  of  f 
zmd  In  f,  as  in  Stewartson  (1974).  The  disturbance  is  now  considered  to  be  proportional 
to 


E  =  exp 


The  disturbance  is  talcen  to  have  slowly  varying  frequency  and  constant  wavenumber, 
so  that  6  2md  are  given  by 


6  —  di  -i-  6^2  T  111  T  "f  •  •  ■ » 


(3.1a) 


and 

fi(t)  =  Qi(t)  -j-  ef]!2(0  "b  liiff^3L(0  “b  f^fl3(0  "b  ■  ■  ■  • 
A  multiple  scales  form  for  the  temporal  derivative  is  introduced, 

d  d  _2  d 
df' 


(3.16) 


(3.2) 


Notice,  that  to  determine  the  effects  of  unsteadiness  the  expansions  need  only  go  as 
far  as  the  fourth  order  terms  not  the  fifth  order  terms  as  in  Smith  (1979  a),  this  occurs 
due  to  the  difference  in  scalings  associated  with  the  temporal  and  sjiatial  variables. 

§3.1  The  Main  Deck 

The  transverse  coordinate  in  this  layer  is  K  =  0(1),  where  Y  —  (~*y.  The  main 
deck  corresponds  to  the  conventional  ^  boundary  layer.  The  l^asic  fiow  in  the  layer 
is  given  by  the  error  function  as  .shown  in  section  2.  It  should  be  noted  that  there  is 
no  normal  component  of  this  velocity,  as  can  be  seen  by  considering  continuity  in  this 
layer  in  the  absence  of  any  disturbance.  It  may  be  shown  that  the  basic  flow  has  the 
properties  that. 


UB{t,Y)^  X{t)Y  +  \3{t)Y^  +  0{Y^)  as  Y 0,  (3.3a) 


UB{t,Y)—*l  as  F  — ♦  oo. 


(3.36) 
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where  it  is  known  that  A(#)  =  At  2,  with  A  =  The  perturbation  to  the  ijasic 

state  in  this  layer  takes  the  form, 

ti  =  [ui  +  (U2  +  f^\n(U3L  +  f^U3  + 

V  =  [fUi  +  f^V2  +  hlfU-JL  +  f^  t>3 

p  =  [fpi  -f-  e^P2  +  hifpsL  +  (^P3  + 

Notice  that  the  logarithmic  terms  (characterised  by  subscript  L),  occur  at  one  order 
higher  than  in  the  corresponding  expansions  of  Smith  (1979  a).  This  is  due  to  the 
form  of  Ufl  as  V'  ^  0  given  by  (3.3  a).  This  form  of  the  disturbance  is  now  substituted 

into  the  governing  equations,  and  nonlinear  terms  are  neglected.  The  operators  A/] 

and  A/2  are  introduced,  and  are  given  by 

do  dxi 

Ml  {u,v)  =  i6iu  +  ,  M2  (u,v)  =  iOiUffU  +  V-:—, 

oY  at 

where  these  are  identical  to  those  used  in  Smith  (1979  a). 

The  equations  in  this  layer  are  conceptually  identical  to  those  of  Smith  (1979  a), 
except  the  order  at  which  the  logarithmic  terms  occur.  In  this  deck  iL'^re  is  no  effect 
from  the  multiple  scales  form  to  the  order  considered,  so  it  is  sufficient  just  to  quote 
the  equations  governing  the  logarithmic  terms, 

dj) 

Ml  {U31,V3l)  +  iO^iUi  =  0,  M2{u31,V3l)  +  UBiOsLUl  =  0,  0=--~. 

oY 

The  boundary  conditions  for  the  disturbances  in  this  layer  are  given  by  matching  with 
the  lower  and  upper  decks,  the  solutions  for  this  layer  are  considered  in  section  3.4. 

§3.2  The  Lower  Deck 

In  this  deck  the  order  one  variable  y  is  introduced,  defined  by  y  =  (~^y.  The  lower 
deck  contains  the  wall  layer  and  critical  layer  of  classical  Orr-Sonimerfeld  theory.  The 
wave  in  this  deck  expands  as 

u  =  [l/i  +  (U2  +  In  fUzi,  +  +  . . .]  E, 

V  =  [,^Vi  +  +  e"  In  (V3L  +  e'Ki  +  ...]E, 

p  =  [f/’i  +  In  (P3L  +  *"^^3  +  . .  •]  £^. 

The  basic  flow  is  given  by  ug  ~  eyA(t)  +  f^y'^As  (t)  Y  ■  ■  ■  as  can  be  seen  from  the 
form  of  (3.3  a).  From  the  normal  momentum  equation  it  may  be  shown  that  the 
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{j  =  1,2,3)  are  independent  of  y,  but  Uj  and  Vj  are  functions  of  y  and  t.  Again 
following  Smith  (1979  a)  the  operators  L\  and  L2  are  introduced,  aiad  are  given  by 


Li  (u, e)  =  i6i u  + 


dv 

w 


L2 {u,v,p^  = i (Ayffi 


—  Qi)  u  +  Au  -h  i6ip  — 


u 

w 


Again  it  is  sufficient  only  to  quote  those  equations  that  diifer  from  those  of  Smith 
(1979  a), 

L,{U,L,V3L)+l03Llh=Q, 


L2  Pii)  —  i^ziUi  -i-  AyW^iL-i  +  i^nPi  —  0. 

^2  (t'^Si  ^3?  -pj)  ~  *  +  ^3^1  )  +  (^2^2  +  )  +  Z  (02P2  +  ^iP\) 

+  iA^SiUiy^  +  Zy^  A^V\  +  =  0. 

at 

The  boundary  conditions  in  this  layer  are  those  of  no-slip  at  the  wall  and  matching 
with  the  main  deck. 

§3.3  The  Upper  deck 

In  this  third  deck  a  coordinate  y  is  introduced,  defined  by  y  =  e^y,  where  y  is  an 
order  one  quzmtity,  and  the  disturbance  takes  the  form 

u  =  [eHi  +  t^U2  +  c^lneusi  -f-  -f- . . .]  E, 

V  =  [evi  -f  -f-  IncvaL  -f  e^tls  -I- . . .]  E, 

p  =  [epi  -I-  c^p2  +  In  cp3L  +  e^P3  +  . . .]  E. 

Here  the  basic  velocity  is  almost  the  uniform  flow  solution  that  is 


Ufl  =  1  +  0(e*). 


In  this  layer  all  the  flow  queintities  are  functions  of  the  slow  time  scale  t  and  y.  The 
operator  D  is  introduced,  defined  by 


D{p) 


d^p 

dy^ 


After  some  algebraic  manipulation  to  remove  the  velocity  components  the  governing 
equations  for  the  pressure  are  obtained  eis, 


E(pi)  =  0,  E(p2)  =  20i^2Pi, 

D  (p3l)  =  2d^03Lpi,  D{p3)  =  29i  (^api  +  92P2)  +  0lp\  ■ 
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The  boundary  conditions  are  those  of  matching  with  the  main  deck  as  y  — >  0  and 
the  requirement  that  the  solution  remains  bounded  as  y  — >  oc.  Having  set  up  these 
equations  they  are  now  solved  in  each  layer  and  the  eigen  i)roblenis  are  derived.  It  is 
worth  stressing  at  this  point  that  there  is  no  effect  in  the  uppt  i  deck  equations  from 
the  unsteadiness  (to  the  order  considered).  The  unsteady  derivatives  are  of  a  higher 
order  unlike  the  spatial  derivative  which  occur  in  the  diffusive  terms,  from 

dx) 


§3.4  Solutions  of  the  main  deck  equations 

The  solutions  of  the  first  order  equations  in  the  main  deck  can  be  shown  to  be 


Vi  =  -tdiAiUB, 


Ui 


Pi  =  -Pi , 


where  Pj  =  Pi  (t)  and  Ai  =  Ai  (t)  are  unknown  slowly  varying  amplitude  functions. 
The  pressure  in  this  deck  will  be  matched  with  the  pressure  in  the  lower  deck.  Solutions 
here  are  the  san\e  as  those  given  in  Smith  (1979  a).  The  integrals  H,  defined  in  Smith 
V.l979  a)  are  functions  of  x  and  Y,  whereas  here  they  are  functions  of  t  and  F,  also 
obviously  ug  is  the  error  function,  rather  than  the  solution  to  the  Blasius  equation. 
The  logarithmic  solutions  are  given  by, 

duB 

J  'W' 

V3L  =  -iOlA-iiUB-,  P3L  =  PsL, 


§3.5  Solutions  in  the  lower  deck 

To  determine  the  behaviour  of  Pi  in  this  deck,  we  differentiate  the  leading  order 
equation  of  continuity  with  respect  to  the  decks  transverse  coordinate,  that  is  y,  and 
eliminate  Vi.  It  is  required  that  the  disturbance  is  bounded  as  y  — >  oo,  thus  the 
solution  for  Pif  is  found  to  be  given  by 


Pi^  =  P,Ai(0, 


where 
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Ai(^)  is  the  Airy  function  and  Bi  is  an  unknown  function  of  t.  Applying  the  no-slip 
condition  at  y  =  0  in  the  leading  order  equation,  it  is  found  that 


(3.4) 


2 

where  subscript  0  corresponds  to  evaluation  at  ^  and  a  prime 

denotes  differentiation  with  respect  to  The  matching  condition  on  the  leading  order 
disturbance  as  y  — >  oo  (therefore  as  |^|  — >  oo),  is  given  by  U\  — >  AAi,  thus  it  is  found 
that 

B\  K  —  AAi , 


where 

OO 

«  =  /  Ai  (t)  di- 

(o 

This  yields  the  first  relation  between  the  pressure  and  the  negative  displacement,  Pi 
and  Ai  respectively,  obtained  by  eliminating  Bi  from  the  equation  (3.4).  Note  that 
these  are  effectively  the  classical  critical  layer  equations,  since,  as  noted  previously  on 
the  lower  branch  of  the  stability  curve  the  critical  layer  coincides  with  the  wall  layer. 
As  can  be  seen  from  Smith  (1979  a)  similar  techniques  may  be  used  to  solve  for  the 
higher  order  equations.  The  same  solutions  as  Smith  (1979  a)  are  obtained  excejjt  the 
unknown  amplitudes  are  now  functions  of  t,  rather  than  x.  The  equation  to  be  solved 
at  the  next  order  is 


d^U2  ,dU2 

ae  ^  di 


—  ® 


dUi 

di 


A"  5  02 


dUi 

di  ' 


which  yields  the  solution 

U2(  =  ^2 Ai  ~^2^i  * BiAi"  +  iociiBi  Ai^ A  ® . 


The  function  Uzl  satisfies 


di^UzL  ,dUzL 
^  di 


—iO, 


A  3  - 


( 


+  iA  + 


A-hzL 


dUi 
di  ’ 


which  gives  the  solution, 

UzLi  =  -S3LAi  +  Ai'"  +  fA-5  Bi  Ai'A-i 
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At  next  order  the  function  satisfies 


de 


j.dU'i  _2  dU2  -rv  A  ^  dU\ 

=  3-^-zQ3A-3 


di 


di 


A<  i  As 


+  -nrA  +6—  _  ,^Vi 


[Qu  aa  a\V. 

TrlfiT’xj  ri))^' 


iX6\ 

j=Q 


t{xe,y 


dUi 


which  yields  the  solution 


dU: 


i=0 


where, 


5Ai  ai  5^Ai 

f  (U  =<*0-^  +  —~^7r  +  02 


+ 


3 

6A3  ^  _  i  /  dAi 


iX9] 


A"3  - 


fld^Ai 

V?  de 

- 1  ^  I 


1  5^Ai  3  ^Ai  ^  „5Ai^ 

2'W"^^~WJ 


\d(  J  J  ^  y  i{X0,y 

where  the  limits  of  the  integrations  are  ^  to  00.  In  this  solution  the  terms  included  in 
the  ^  summation  include  those  of  Smith  (1979  a)  and  the  unsteady  terms.  The  terms 
included  in  the  function  /((f)  eire  those  eissociated  with  A3,  that  is  the  second  order 
mean  flow  quantity. 

The  /3's  are  given  by, 

0  2  A 

00  =  B3,  01  =  aiiB2  +  0112B1  —  -^Biaii  +  A“3  ’ 


As  dBi  2d  p 


03 


=  i  (7S.  +  ^  («2B2  +  «3B,  )  -  , 


*  =  -4^' 


1  /02  \ 

^4  =  ^  I  +o;n«2  I  , 

where  7  has  been  introduced,  and  is  given  by 

^  9iX\ni  X)' 


00  = 


02(^2 

W’ 


12 


The  terms  ao,  oi ,  02  and  03  are  given  by 


A3  .  _2  A3  _1  3Qi 

a„  =  -A  3,  =_A  3—,  a. 


_  3A3  /  *  1  /  ^^1  \ 

~  ^  VA^i  j  ’  ~  \  ^  VA^i  j 


For  brevity  ai  j  and  02  have  been  introduced,  and  axe  given  by 

aij  =  -iA~5  and  a2  =  -^Bi(0- 

Now  applying  no-slip  conditions  at  the  wall  it  is  found  that 

i{e,P2  +  ^2^1)  =  +  A3;^B,Air  -f-  onBiAi'', 

OU\ 

ii^iPzi+BziPi)  =  AsBaiAi;  -f  A^^BiAir  +  *A-t  BjAi", 

i(0iP3+62P2  +  e^Pi)  =  Ai 

i=0  ^ 

Matching  with  the  main  deck  yields 

B2K  —  0'2Ai()  —  cvii^iAio  =  A  ^^2  —  P\I  —  , 

Bslk  -f  Ai'o'  -  iA't  -  03l)  BiAio  =  A  (a,l  -  PJ  -  , 

fi 

_2  ^— ^  ^  A  ^  f  r  \  (  A  rt  f  f  ^2  i  ^2  A  ^3  A  \ 

B,k-A-3  ^ftA.0  -A3  j  fd(  =  X  {A,  -  P.,1  -  -j—I  -  +  ^A,  -  j  . 

(3.6) 

The  integrals  occurring  in  the  expressions  (3.6)  are  given  by 

i  =  ^_ul\Y,)dYu  and  I  =  {Y,)dY,, 

where  A  is  arbitrary  and  non-zero.  The  horizontal  bar  denotes  that  only  the  (Hadamaxd) 
finite  part  is  to  be  retained.  We  note  that  for  small  77  the  error  function  can  be  ap¬ 
proximated  by,  A  (t/  —  r;*/3  -b  and  hence  I  is  given  by 


T  ^  1 

1  ^  ~  T*  — 

A2  A 


85A  -  12A^ 


A  26  (30  -  10A2  +  dX*^ 


1  r  43 
26  7  30- 


435  -  907/2 


10t;2  -|-  3t;^ 


dr]  , 
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with  a  similar  expression  for  I.  In  the  numerical  calculations  the  quantity  H2oo  —  H\oo  — 
i  is  found  to  be  — 0.29512(2t5  this  is  required  in  the  second  order  eigenrelation  given 
by  (3.8  b). 

§3.6  The  upper  deck  solutions  and  eigenrelations 

The  pressure  equation  in  the  upper  deck  yields  the  bounded  solution 

where  P\  —  P\  (<)  is  an  unknown  amplitude  function.  Now  matching  with  the  main 
deck  as  y  — >  0  yields 


Pi— Pi,  6\P\  —  9\A\  id\V\oo)  1 


where  uioo  denotes  the  values  of  v\  (the  main  deck  normal  velocity)  as  the  decks 
transverse  coordinate  (T )  tends  to  infinity.  This  provides  the  leading  order  eigenvalue 
problem 

=  XA^Ai'o.  (3.7) 

This  relationship  determines  the  wavenumber  in  terms  of  the  frequency  and  the  skin 
friction  A.  The  higher  order  pressure  functions  may  be  determined,  and  are  given  by 

P2  =  (A-^2Ply)e-‘'>^ 

P3L  =  {PzL  —  ^SiPiy) 

P3  =  [pz  -  {BzPi  +  92P2)  y  +  ^^2Piy') 

To  match  with  the  main  deck,  the  unknown  functions  must  satisfy 

A=P2-A,02^1oo, 

9i 

i9\P\H2oo  d"  2if2i.<4  —  i9iA2  =  iA\92  —  iPj  —  i— Pi, 

9i 


P3L  —  PzLy 

-9\A^i  =  i9^iA]  —  iP^L  - 


iBziPi 
9,  ’ 


P3  —  P3  +  9\P]Hzoo  —  9i  {9\A2  +  ^2-^1) -ffioo  +  2Q,\9\A\Hiooi 
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ii02Pi  +  ^1^2) ^^200  -  tOlAiH^oc  +  iniPiHeoo  +  (a2  - 

+  —  i^iAa  =  —  ((fii  —  ^2)^200  +  (^^2  —  ^3)^100)  —  iPz  —  i  \  ^ 

"1  t’l 


^2-^2 


The  higher  order  eigenrelations  may  now  be  obtained, 

DBo-C 


AipP 

kA5 


D631  —  ilsLdi 

^  AioD 


/cAa 


(^200  —  Hioo 

(3.8a) 

AioD 

kaI  “  A4  ’ 

(3.86) 

i(c.t.a  +  C.t.Xa  +  n.t.)  , 

(3.8c) 

where  D  and  D  are  given  by, 

D  =  l+^, 

Ain 


In  (3.8)  c.ta  denotes  the  terms  that  arise  in  conventional  analysis;  these  are  exactly 
the  same  as  those  occurring  in  the  third  order  terms  of  Smith  (1979  a),  as  the  non¬ 
parallelism  has  no  effect  at  this  order.  These  terms  can  be  directly  determined  from 
(3.15  b)  of  Smith  (1979  a).  The  terms  c.f.Aa  axe  those  associated  with  the  second 
order  basic  flow  terms  and  n.t.  are  the  terms  that  arise  due  to  the  f-derivatives  (that 
is  unsteadiness  of  the  modes).  Discussion  of  the  solution  of  these  eigenrelations  is 
included  in  section  5. 


§4  Nonlinear  theory  for  larger  disturbances 

We  start  by  setting  up  the  nonlinear  problem  for  an  unsteady  triple  deck.  As  in 
the  linear  problem  the  scalings  x  =  and  t  =  e^T  are  introduced.  In  the  Rayleigh 
layer  under  consideration  at  the  time  station  t,  the  flow  field  in  the  absence  of  any 
disturbance  is  given  by 

V  =UBit,Y)  +  ...  , 

where  UB{t,Y)  is  given  by  the  error  function.  It  is  also  worth  re-iterating  at  this 
point  that  as  F  — >  0, 

ub(<,F)-^  A  (t)y-F  A3  (t)Y^  +  ---, 
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where  X{t)  =  (7r<)~^. 

The  properties  of  the  entire  problem  may  now  be  determined  by  examining  the 
lower  deck,  the  same  transverse  coordinate  is  used  as  in  section  3.2,  that  is  y  =  f^y, 
where  y  =  0(1).  In  the  lower  deck  the  flow  is  taken  to  be  of  the  form, 

u  =  {eU  (X,  y,  Z,  T;  t) ,  (X,  y,  Z,  T;  t) ,  eW  (X,  y,  Z,T;t)  P  (X,  Z,  T;  t))  . 


This  flow  includes  both  the  basic  flow  and  any  disturbance  quantity.  The  fact  that 
the  pressure  P  is  independent  of  y  is  given  by  the  wall  normal  momentum  equation. 
The  governing  equations  in  the  deck  axe  now  found  to  be 


dU 

rrdU 

.rdU 

dP 

d^U 

dT 

dy 

4-  W -  =  - 

dZ 

~  dX  ^ 

dr' 

dW 

.rdW 

dP 

d^W 

dT 

+  F  — 
dy 

+  ”^82  = 

dZ 

dr 

dU  dV  dW  ^ 

dX  dy  dZ 

with  boundary  conditions 


(4.1a) 


(4.16) 


(4.1c) 


U  =  v  =  W  =  0  at  y  =  0. 


(4.1d) 


eind. 


U  ^  X{y  +  A{X,Z,T))  as  y  ^  oo. 


(4.1e) 


Here  it  is  assumed  that  the  spanwise  component  of  the  disturbance  is  confined  to  the 
lower  deck.  From  the  upper  deck  equations  we  see  that 

do  d^A 

=  0,  with  p  =  P  (X,  Z,  T)  and  —  =  at  y  =  0 

where  y  is  the  transverse  upper  deck  coordinate,  and  p  is  a  function  of  X,  y,  Z  and 
r,  and  is  bounded  as  y  — ^  oo.  Now  consider  a  small  perturbation  to  the  basic  state 
(U  =  Ay,  V  —  W  =  P  =  A  =  Qi)  'm  the  lower  deck.  The  perturbation  to  this  system  is 
taken  to  be  of  order  1.  In  Hall  h  Smith  (1984)  initially  a  bi-modal  analysis  was 
considered  and  subsequently  a  multiple  mode  analysis,  here  a  single  mode  is  considered. 
It  is  assumed  that  h  is  larger  than  any  positive  power  of  e  so  the  governing  equations 
can  be  assumed  to  be  given  by  (4.1)  for  all  orders  of  h.  In  section  6,  the  problem 
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involved  in  malcing  /i  ~  e"  is  discussed,  where  unsteady  effects  occur  at  the  same  order 
as  nonlinear  effects.  The  flow  is  now  perturbed  and  taken  to  be  given  by 

U  =  Xy  +  hUi+  h‘^U2  +  h^U3+...  , 

V  =  hV,  +  h'^V2  +  h^V:i  +  . . .  , 

W  =  hWi  +  h‘^W2  +  +  . . .  , 

P=hP,+h‘^P2  +  h^P3+...  , 

A  =  hAi  +  A2  +  h^A^  +  •  •  •  , 

p  =  hpi  +  h'^p2  +  h^pz  +  ...  , 

The  disturbance  is  considered  to  be  proportional  to  E  and  its  integer  multiples,  where 
E  is  given  by 

E  =  exp  (i  {9X  +  7Z  -  ilT)) , 

7  and  Q,  are  real  constants  for  neutral  stability,  and  Z  represents  the  spanwise 

_  3 

variation  on  the  0{Re  ®)  scale.  For  any  spanwise  wavenumber  7  a  pair  6  and  may 
be  determined  using  the  first  order  eigenrelation  given  by  linear  theory.  At  this  point 
a  new  time  coordinate  is  introduced,  namely  T  defined  by 

T  =  — T  =  h~^T 
111 

This  necessitates  the  inclusion  of  further  terms  in  the  skin  friction  in  the  neighbourhood 
of  the  neutral  time  <oi  as  mentioned  in  Smith  (1980);  the  correction  to  Smith  (1979 
b),  so 

t  —  to  A  h^i,  A  =  A  +  \2  "I"  ft  =  fti  -H  /i^ft2  0(/i^), 

where  A2  =  tidX/dt)\t=to-  If  is  necessary  to  include  the  further  terms  in  the  ft  form, 
to  satisfy  the  third  order  eigenrelation. 

Now  substitute  the  disturbance  into  the  governing  equations  including  the  neces¬ 
sary  multiple  scales  approach  for  temporal  derivatives.  The  equation  of  continuity  at 
successive  powers  of  h  becomes, 

Li  (l/i ,  Fi ,  tFi )  =  0,  Li  (C/2,  V2,W2)  =  0,  L,  (C/3,  V3 ,  =  0. 

The  streamwise  momentum  equation  at  successive  orders  yields, 

L2{UuVi,P,)  =  0, 
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,  „du,  du,\ 

LAU,.Vi,P^)  =  -[^U,g^+V,gy+W,-^j, 

i2  (V,,V,,P,)  -  -  -^  -  -  r.A; 

dVi  dU,  dU2  dUi 

-'''W-^^W-'^'lz-''^-dz- 

The  spanwise  momentum  equation  become, 

L^iWuPy)  =  0, 

r  rr.r  r.  ^  rr  ^Wi  51^1  dWi 

Li(W2,P2)-  Ui  Vi  ^  dZ  ' 

dH%  dWt  dW2  ^W^ 

*  dY  ^  dY  ^  dZ  ^  dZ  ■ 

The  operators  L\,  L2  and  L3  are  defined  by, 

r  /rrTrr.rx  dV  dW 

Li  (U,v,w)  -  —  +  ^ 

T  /rr  Tr  A  VT.  ^P 

2(U,  ,P)-^  +  dX  dY'^' 

Lii  -  or  QX  dZ  dY^‘ 

The  boundary  conditions  at  successive  orders  are  given  by 


Ui  =  Vi  =  Wi  =  0  at  y  =  0  (i  =  1,2,3), 


Ui  —y  XAi  (X,Z,T)  as  y —y  00, 

U2  —y  XA2  {X,Z,T)  as  y  —y  00, 

U3  ^XA3(X,Z,T)  +  X2Ai(X,Z,T)  as  y -y  00, 

It  is  worth  noting  another  significant  discrepancy  between  this  work  and  that  con¬ 
sidering  a  spatially  growing  boundary  layer,  that  is  the  slow  time  derivative  does  not 
come  into  the  upper  deck  governing  equations,  as  was  noted  in  the  linear  section.  Now 
consider  Ui  to  be  wavelike  proportional  to  E  and  its  conjugate  denoted  by  E.  The 
first  order  disturbance  field  is  now  given  by 

Ui=Ui(y,f')E  +  c.c.,  Vi=Vi(y,f)E  +  c.c.,  Wi  =Wi  (y,T)  E  +  c.c., 
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P\  =  P\  E -k- C.C.,  Pi  =  pu  (y,f^  E  +  C.C.,  Ai  =  Ai  (t^  E  +  C.C., 

where  c.c.  denotes  complex  conjugate  and  the  disturbance  quantities  with  tildes  axe 
independent  of  X,  Z,  and  T.  By  combining  the  first  order  equations  the  two  governing 
equations  axe  obtained 


de  aS  «A  ’ 


d^u,  ^dii, 


where 


?  =  as(s-|),  a  =  .ao. 


(4.2a, 6) 


(4.2c) 


At  this  point  it  is  advantageous  to  introduce  a  further  disturbance  quzmtity  that  is 
Xij-  The  equations  (4.2  a)  and  (4.2  b)  may  be  manipulated  and  it  is  found  that  yn 
satisfies 


where 


^Xii  ^ 


=  iQUi  + 


The  solution  of  (4.3  a)  together  with  a  bounded  condition  at  infinity  yields, 

Xu=Bi  (t)  Ai(0. 

The  three  dimensional  eigenrelation  may  now  be  obtained,  namely 

A2Ai'(eo)  =  (ieA)3  (02^72) 

as  given  in  Hall  &  Smith  (1984),  where 


^0  =  -  ■ 


5  r 

—  amd  K  =  /  Ai  {q)dq. 

J 


(4.3a) 


(4.36) 


(4.3c) 


(4.4a) 


(4.46,  c) 


The  neutreJ  solution  of  the  eigenvalue  problem  is  given  by 

Co  =  -2.298i5. 

In  this  relationship  6  is  real  and  (4.4  a)  yields 

^5  (02  ^72)5  =  I.OOIA^ 


=  2.299^5  A5. 


Thus  for  a  particular  spanwise  wavenumber  7,  a  pair  B  and  flj  maybe  obtained  such 
that  neutrality  is  assured,  a  plot  of  the  variation  with  7  of  0  and  fli  is  given  in  figure 
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Figure  3:  Variation  of  6  and  Qi  with  7. 

The  solution  of  the  second  order  problem  is  now  considered,  the  disturbance 
quantities  now  contain  second  harmonic  terms,  the  disturbance  is  taken  to  be  of  the 
form 

U2  =  EU21  +  +  E^U22  +  +  U20, 

with  similar  expansions  for  Vji  ^2,  P2,  P2  and  ^2-  The  terms  proportional  to  are 
considered  and  it  is  found  that  ^^22?  V22,  IV22, -^22^  satisfies  the  system 

dV22  .  ~ 

2idU22  H — — I"  217^1^22  =  0, 
oy 

^2  rV  ^TT 

—2iQU22  +  \y2i6U22  +  -^1^22  +  2idp22 - =  —i8Uf  —  V}  - ijW\U\, 

oy^  oy 

-2iQW22  +  \y2ieW22  +  2^7^22  -  ijWl 

along  with  boundary  conditions 

U22  —  ^22  —  IV22  =0  at  y  =  0, 

t/22  AA22  as  y  — >  00. 
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Again  at  this  point  it  is  advantageous  to  re-introduce  x,j,  this  time  defined  by 


X22  —  2iBU22  +  2i')W22-i 


so  the  system  becomes, 

^^X22 

dC 


-2^ 


dx 
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dxu  ,  5^X11 

-\XU-g^+V,-^ 


The  solution  of  which  necessitates  the  introduction  of  given  hy  ^  =  2"5^.  Note  that 
it  is  not  necessary  to  solve  the  equation  to  determine  ^C/21 ,  i^2i ,  W^2i  1  ,  as  it  can  be 

seen  that  the  governing  equations  for  these  quantities  are  the  homogeneous  equations 
(4.2);  the  solution  is  a  lineax  multiple  of  the  first  order  solution.  Now  finally  at  this 
order  the  mean  adjustment  components  are  determined  by 


5^20  _  „ 
dy 


t7(0 


dUi 


dp  ' '  dy  dy 


d^W: 
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dp 

with  boundary  conditions 


, 


=  -ie 


U20  —  V20  =  W20  ~  0 


at 


y  =  0. 


dy 


By  the  usual  substitution,  that  is  X20  =  iBU2o  +  *714^20,  it  is  shown  that 

dV2o 


di 


=  0, 


^(xnXu)  ■ 

Finally  the  third  order  equations  are  considered  and  ageun  by  elementary  analysis  it  is 
known  that  the  third  order  quantities  contain  terms  proportional  to  E^,  E^,  E^ ,  E°, 
E^ ,  E^  and  E^,  however,  determination  of  the  functions  associated  with  E  is  sufficient 
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to  give  the  required  amplitude  equation.  The  functions  V31,  If  31 ,  are  give 

by  the  solution  of  the  system 


1  du,  ,_,7 


+  um) 


“aT  '“ST  “af 


+i7  (2C22IC;'’  -  W22UI"'  +  W'loW'i))  , 

+  ^  (21^22^1''^  -  +  ^20 1^1) 

A  a  \  / 


+  ^  (21^22^1''^  -  ^22f^;‘^^  +  ^20 1^1) 

J_  { i>(c)^  (7.  ^  +  Vio^l 

AM  ‘  di  ^  di  ^  di  di 


(^22 


Wi^'*  +  W20W,)  , 


with  boundary  conditions 

U31  =  V31  =  W31  =0  at  y  =  0, 

1^31  ^  AA31  +  A2A11  as  y -e  00. 

We  meike  the  substitution  X3i  =  i^^si  +  *7^31  and  it  is  found  that 


^^X31  _  ^^X3i  ^  .  -2C)^XU 
nc3  ^  nc 


-I  ^Xii 


d^dT  XV  XdJ 

As  _  9x22  AaTr^^X20  As^-^a^x 


+  1  -ir-Xii 


^Ur^X20  IS--,U-X22  .i,,  C7-X1] 


,  A  4  ^Xn  A  A  ^Xn  \ 

+As-^^X22  -  As-^-X2o  1 


We  require  a  solvability  condition  for  this  equation  (with  the  upper  deck  terms),  which 
yields 

^  =  a^Sl  +  ^2Bl|5l|^  (4.5) 
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»vhere  B\  is  the  lower  deck  amplitude  function,  and  where  is  an  integral  gi'.t'ii 
in  the  appendix  (equation  (A.l)).  Note  that  <7i  may  be  thought  of  as  a  frequency 
perturbation,  and  if  a\  =  ifi,  we  may  determine  (l  from  the  eigenrelation  (4.4  a), 
along  with  (4.4  b,c).  If  J7  is  replaced  and  A  bj'  A  +  X  {(l  —  and  A  =  ). 

we  obtain  the  relationship, 

4-  (0'^  +7^)' 

(B-  +  -  As^o) 

Note  that  it  is  not  necessary  to  calcidate  the  vahu'  of  a-/  as  it  has  the  same  value' 
as  in  Hall  &  Smith  (1984),  where  it  is  bt'e'ii  assuinerl  that  A  is  unity.  The  velocity 
field  has  been  norniali.sed  so  that  ^^|v  :=o  =  1-  These  results  are  for  a  single  modt', 
whereas  Hall  &:  Smith  (1984)  considered  two  modes,  note  that  the  value  of  <iir  (vdiich 
corresponds  to  here)  is  independent  of  the  secoml  spanwise  wavenumljer  thus. 


r 

■a 

3.9 

3.6 

3.3 

"1 

3.0 

1 

2.4 

2.1 

1.8 

1.5 

1.2 

H 

-0.0039 

-0.0049 

-0.0081 

-0.0147 

-0.0214 

-n.0314 

-0.0515 

-0.0919 

-0.2250 

Table  1,  7  against  a2r,  as  in  Hall  &  Smith  (1984). 


Note  that  Q  =  q]  (i  —  t„)  and  A  =  X'  (to)  (i  —  t ,  hence  the  relationship  (4.6) 
becomes 

A'(fo)  (i.46912A5) 

A  (to)  (2.29726>5  (B^  +  7^)^  +  5.266128A5)  ’ 

where  a  superscript  prime  represents  differentiation  with  respect  to  f.  It  is  noted  that 

A'  (To)  <  0,  thus  Qi  >  0  after  the  neutral  time  i  =  tn,  as  is  expected  from  linear  theory. 

Note  that  by  considering  smaller  disturbances  unsteady  effects  may  be  incurred. 

It  is  assumed  that  h  is  small,  but  not  so  small  that  the  leading  order  effects  described 

in  this  section  are  changed.  In  fact  the  calculations  included  in  this  section  are  appro- 

_ l_  _ L  , 

priate  to  1  ^  /i  i2e  •  The  case  where  /j  ~  '®  is  now  considered,  note  that  it 

is  simila’’  to  the  case  discussed  in  Hall  &  Smith  (1984),  for  the  spatial  problem.  It  is 
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taken  that  h  =  HRe^ ,  with  H  of  order  one,  and  corresponding  to  a  disturbance  size 
—  * 

of  ^  =  HRe  ^ ,  and  it  is  found  that  the  nonlinear  amplitude  equation  is  given  by, 

^  =  iqi  (-U  +  i  +  TH-*)  B,  +  (725, 15, 1^  (4.7) 

dT  '  ' 

The  coefficients  in  this  equations  remain  the  same  as  in  (4.5)  so  no  further  numerical 
work  is  needed,  the  discussion  of  the  solutions  of  this  equation  are  included  in  section 
6. 


§5  The  neutral  stability  characteristics 

§5.1  Quasi  steady  flow  stability  in  the  triple  deck  (linear  case) 

Here  we  find  the  wave  characteristics  that  can  be  determined  from  the  eigenre- 
lations  (3.7  h  8).  It  should  be  stressed  at  this  point  that  triple  deck  theory  is  only 
applicable  for  Re  »  1,  however  we  can  compute  the  neutral  wavenumber  6n  and 
compare  it  with  the  Orr-Sommerfeld  predictions.  The  eigenrelations  (3.7,8a  &:  b)  are 
unaffected  by  the  boundary  layer’s  growth,  so  the  condition  of  neutrality  is  realised 
by  ensuring  that  fi,,  and  are  real. 

To  ensure  that  Q.\  and  6\  are  real  we  calculate  the  value  of  (defined  in  section 
3.5),  required  to  make  the  quantity  Aio/(iiK)(=  c,  =  1.00065)  real.  This  occurs  when 
^0  ~  — 2.298z5  =  — which  yields 


3  ^ 


=  Af  eocf  =  Af  n, ,  =  Ai^ c’  =  A’ By . 


At  this  point  we  recall  that  but  $1  is  a  constant  so  we  can  think  of  this 

expression  as  defining  the  neutral  time, 


After  this  time  modes  will  exhibit  temporal  amplification,  and  the  weakly  nonlineeir 
anEilysis  described  herein  will  come  into  effect  eis  the  mode  attains  a  finite  amplitude. 
It  appears  that  the  higher  wavenumber  modes  will  start  to  grow  earliest,  of  course  this 
can  be  seen  from  figure  1.  These  results  may  also  be  interpreted  in  terms  of  the  local 
boundary  layer  Reynolds  number,  defined  by  Rs*  =  2/y/7r{tRe)^ ,  where  2y/i/y/TT  is 
the  displacement  thickness.  As  noted  in  section  2,  the  Orr-Sommerfeld  calculations 
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are  affected  for  t  =  1/4  so  that  the  value  for  the  displacement  thickness  is  \\^ 

can  express  the  neutrally  stable  wavenumber  in  terms  of  this  Ren  nolds  number  as 


n 


25c? 

5 

TT  4 


We  now  proceed  to  the  second  order  eigenrelations,  and  note  that 


Cl 


2^1  Aio 
2  — 


We  require  that  0,2  and  $2  are  real,  and  thus  we  find  that 


62  =  1.55721t~5. 


In  order  to  determine  the  logarithmic  wave  characteristics  we  note  that  A3  =  —  l/(12\/7rt^), 
and  we  find  that 

$31  =  -0.56748f5. 

Subsequent  terms  in  the  expansion  are  effected  by  the  growth  of  the  layer,  hence  we 
defer  their  discussion  to  the  next  section.  We  may  combine  the  expressions  for  B\ ,  62 
and  $31  in  order  to  predict  the  neutral  wavenumber  in  terms  of  the  local  boundary 
layer  Reynolds  number.  We  find  that 

=  1.0789/27^  (1  +  3.37079^7/  -  5. 70354777^  lni7«. 

this  value  is  compared  with  the  neutral  value  found  from  Orr-Sommerfeld  theory,  the 
results  are  displayed  in  figure  4.  In  order  to  predict  the  upper  branch  we  would  need 
to  consider  a  more  complicated  structure  in  which  the  critical  layer  is  distinct  from 
the  viscous  wall  layer. 
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Figure  4:  Comparison  of  neutral  wavelength,  Orr-Sommerfeld  and  triple  deck 


The  agreement  is  reasonable,  but  it  should  be  noted  that  analysis  relies  on  different 
powers  of  the  Reynolds  numbers  being  distinct,  however  at  the  values  of  Rs>  displayed 
in  figure  4,  Rg.^  IniZ^.  and  Rg»*  are  not  totally  distinct.  However,  the  triple  deck 
theory  does  appear  to  capture  the  lower  branches  characteristics. 

§5.2  Unsteady  effects  on  triple  deck  stability 

In  this  section  we  discuss  the  effects  that  the  boundary  layer’s  growth  has  on 
the  stability  of  the  Tollmien-Schlichting  waves.  As  heis  been  mentioned  several  times 
earlier,  the  evolution  of  the  layer  changes  the  characteristics  at  one  order  higher  than 
in  the  spatial  case  of  Smith  (1979  a),  thus  it  is  perhaps  even  more  crucial  in  these 
temporal  cases  that  the  layers  evolution  is  included  in  the  analysis.  This  is  to  be 
expected,  since  during  one  spatial  wavelength  of  the  wave  the  boundary  does  not 
evolve  (spatially)  as  much  as  it  would  during  one  period  (temporally).  The  leading 
order  effect  of  the  boundary  layer’s  growth  occurs  in  the  eigenrelation  (3.8  c),  which 
determines  fla  and  $3. 

We  consider  the  total  temporal  derivative  of  a  flow  quantity  $,  and  define  its 
growth  rate  as, 

1  5$ 

$  at  ■ 
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We  require  that  the  quantity  cr  is  purely  imaginciry,  which  for  a  problem  which  does  not 
exhibit  any  temporal  changes  translates  to  ensuring  that  the  fl’s  are  real.  However, 
we  know  that  dj  dt  +  so  that 

£i 

$ 

Hence  in  order  for  a  to  remain  totally  imaginary  (to  the  order  considered),  we  require 

Imag(fi3)  =  —Real 

Using  this  condition  together  with  eigenrelation  (3.8  c),  we  can  determine  9^  and 
[Real(f23)  -  Imag(^)]. 

It  is  not  clear  a.s  to  which  quantity  should  be  used  in  order  to  represent  the 
disturbance’s  growth.  In  Smith  (1979  a)  the  v'all  pressure  was  used,  and  the  results 
which  were  obtained  allowed  enhanced  agreement  with  the  experiments  of  Ross  et 
al  (1970).  Smith  noted  that  this  implied  that  the  boundary  layer’s  growth  had  a 
stabilising  effect,  however  if  other  quantities  were  chosen  to  represent  the  disturbance, 
the  conduaions  could  be  the  opposite.  In  Gaster  (1974)  and  Eagles  &  Weissman 
(1975),  it  was  noted  that  the  difficulty  is  in  defining  the  instability  in  some  absolute 
sense,  however  this  method  should  allow  good  agreement  with  a  particular  experiment 
(in  which  the  disturbance  is  mea.sured  using  a  certain  quantity).  We  are  not  aware  of 
any  experimental  data  concerning  the  susceptibility  of  a  Rayleigh  layer  to  Tollmien- 
Schlichting  waves.  If,  for  example  we  take  the  wall  layer  normal  velocity  (Vj)  to  be 
representative  of  the  disturbance,  then  the  layers  growth  is  seen  to  have  a  destablising 
effect. 

The  reason  for  the  enhanced  agreement  with  experiment  seen  in  Smith  (1979  a), 

—  3 

is  likely  to  have  its  origin  in  the  size  of  the  correction,  it  is  now  of  order  Re  *  which 

_  I 

is  larger  than  that  of  Bouthier  (1973)  and  Caster’s  (1974)  (order  Re^)  successive 

_  I 

approximation  technique.  As  mentioned  previously  our  correction  is  of  order  Re  “ 
which  represents  an  even  larger  change,  so  we  expect  even  better  agreement  between 
the  theory  and  experimental  results. 


<T  =  e  ^  — zfii  —  fifl2  —  In  (ilsi  —  e  (  — 
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§6  Discussion  of  the  weakly  nonlinear  stability. 

The  relationship  between  the  equilibritun  amplitude  solutions  of  (4.5)  and  those 

of  the  unsteady  flow  amplitude  solutions  of  (4.7)  is  now  considered.  The  latter  as 

_ 1_ 

previously  stated  corresponds  to  a  disturbance  of  size  /i  ~  iZe  The  former  is  for 

—  ‘ 

larger  disturbances  where  h  ^  Re^ .  Thus  we  are  required  to  solve, 

^  =  iClB  +  a2B\B\\  (6.1) 

at 

and, 

^  +  H-*f)  B  +  a2B\B\\ 

where  B  —  Bi  is  the  lower  deck  amplitude  function  as  in  (4.3  c) 

(6.1)  has  the  stable  equilibrium  solution, 

|B|^  =  — ,  (6.3) 

<72r 

which  bifurcates  supercritically  at  Qj  =  0  from  the  zero  solution  for  increasing  Reynolds 

number,  where  a2r  is  always  less  than  zero  as  shown  in  table  1.  Increasing  the  Reynolds 

number  corresponds  to  progressing  through  time,  thus  increasing  to .  The  zero  solution 

is  unstable  beyond  ft  equals  zero,  so  any  perturbation  to  this  tends  to  the  solution  (6.3). 

—  1 

Now  note  that  sufficiently  close  to  the  neutral  time  t„  (specifically  within  h  ^  Re  ^ 
the  appropriate  equation  is  (6.2).  Drawing  on  the  conclusions  of  Hall  &  Smith  (1984) 
we  let  B  =  pexp{i<f>),  where  p  and  <f>  are  real  functions  of  T.  It  was  shown  that  the 
large  T  behaviour  of  p  is  given  by, 

2  ^  \ 

p  ~  -  I  as  i  — >  oo  I  . 

<72r  ^ 

It  may  be  seen  from  the  numerical  results  of  section  5  that  qi  and  o’2r  are  both  always 

negative.  This  equation  implies  that  the  large  T  behaviour  may  be  inferred  from  (6.2) 

neglecting  the  dB/dt  term.  This  equilibrium  solution  corresponds  to  letting  ft  0  in 

equation  (6.1),  wh:'i-ii  corresponds  to  the  disturbance  size  being  decreased.  Thus  the 

interval  over  which  unsteady  terms  induce  a  contribution  to  the  eunplitude  equation 

_  1 

is  crucial.  In  this  interval  a  small  disturbance  {h  ~  i2e  ^),  develops  smoothly  into  a 
finite  amplitude  Tollmien-Schlichting  wave.  Only  ‘later’  (that  is  for  relatively  large 
time)  does  the  amplitude  tend  towards  the  amplitude  predicted  by  (6.1). 


§7  Conclusions 

We  have  shown  that  in  theory  a  Rayleigh  layer  is  unstable  to  Tollmien- Schlichting 
waves.  An  Orr-Sommerfeld  calculation  is  performed  and  it  is  found  that  beyond  a 
certain  Reynolds  number,  disturbances  with  fixed  spatial  wavenumber  comparable  to 
the  displacement  thickness  are  temporally  unstable.  As  the  Reynolds  number  increases 
further  the  structure  of  the  neutral  modes  can  be  described  by  a  triple  deck  method, 
and  eigenrelations  are  calculated  and  used  to  determine  the  wave  characteristics.  A 
multiple  scales  technique  is  used  to  determine  the  effect  of  the  growth  of  the  layer  on 
the  v'aves’  stabilities  In  summary  it  is  found  that  this  depends  on  which  quantity 
is  chosen  to  represent  the  waves  intensity.  A  Stuart-Watson  amplitude  equation  was 
derived  in  the  neighbourhood  of  the  neutral  time,  and  the  waves  were  found  to  bifurcate 
supercritically  from  the  zero  solution.  An  amplitude  is  found  that  represents  a  crucial 
stage  in  the  development  of  the  wave,  this  pahse  allows  the  wave  to  evolve  smoothly 
from  its  linear  stage  to  its  equilibrium  amplitude. 

Although  this  work  in  this  paper  represents  fairly  minor  modifications  to  various 
earlier  papers;  the  methods  are  applied  to  an  entirely  different  physical  problem.  It 
will  be  of  interest  to  determine  whether  the  upper  branch  modes  can  also  be  compared 
with  the  Orr-Sommerfeld  calculation.  The  calculation  will  be  solved  using  the  methods 
discussed  in  Bodonyi  &  Smith  (1981),  wherein  the  effect  of  non-parallelism  on  upper 
branch  Tollmien- Schlichting  waves  in  a  Blasius  boundary  layer  is  discussed.  The  upiJer 
branch  problem  is  far  more  complex,  since  the  critical  layer  is  distinct  from  the  wall 
layer.  Within  the  context  of  this  spatially  homogeneous  problem  it  will  be  interesting 
to  find  what  acoustic  radiation  may  be  generated  by  these  waves,  and  whether  the 
waves  on  the  upper  branch  or  the  lower  branch  of  the  neutral  curve  produce  the  most 
sound.  The  noise  can  be  calculated  using  methods  similar  to  those  discussed  in  Tam 
&  Morris  (1980). 
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Appendix 

The  double  integral  a 2  involved  in  the  weakly  nonlinear  calculation  occurring  in  equa¬ 
tion  (4.5),  is  given  by 


1  I  (02  ^-.2)  2  .  Q  I 

<^2  = - ^ ^ /  LGiq)dq  +  —LG{(o)  , 

^(^2+  2)5  Ai(^o)-Ai''(eo)  I  7  d(  j 


(A.l) 


where  the  function  G(i^)  is  given  by, 


-  .  •(  (if-  r■^  if  f- 

G(0  =  A  3  A  3  <T(a)  I  -  y  Ai  (9)  d^Ai  (^0)  +  2  y  j  k\{q)dqdq — 


J  ki{q)dqAi{i)  +  j  J  Ai{q)dqdq^~ 


r  '  ^ 

^  J  Ai(q)dqL2{FiO)+^  J  J  Ai(9)^i2(F(0) 


^  J  L2iFiq))dqAiiO+^  I  J  L2{F{q))dq^ 


-  A-3 


'  ^  2 

J  J  Ai{0dqd(^  +  Ai{0H{0  , 


2I  (02  +^2)i  ^ 

^(a)  =  - - 777 — 

4(^2  +72)2  ±  J  Ai  {q)dq  +  2i  AsAi' 

?0 

and  the  form  of  F{()  included  in  the  integral  a^a)  is  given  by. 


e(e)  =  ^  /  Ai(q)dqAi(0- J  J  Ai{Od(d(^  , 

\?0  €o  ?0  / 
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the  function  H(()  is  given  by, 


^  «  9 

=  13  J  j  ki 


{q)dqdi 


^0  ia 


dk\ 

'W 


{q)dqdi 


dki 

'W 


(o  (o 


+  ^(Ai(0Ai(0)V 


where  the  operators  L  and  Lj  are  defined  as, 

g(2 

d^6  dd> 

The  terms  c.t.x^  occurring  in  equation  (3.8),  that  is  those  due  to  the  second  order 
basic  flow  terms,  are  given  by 


c.t.xa  = 


where  /(^)  as  given  in  section  3.5. 

The  terms  due  to  non-conventionaJ  analysis  in  equation  (3.8),  are 


n.t. 


=  ir-8rr»-A5r-‘'^'»  )  +  n7'ar(-^  "a.„— -n.A.„  j 


+ 


it  ,AUi;'  + Air  A-^)  +  ;VA»Air 


{A.2) 
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